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We reveal extraordinary electromagnetic properties for a general class of uniaxially polarizable 
media. Depending on parameters, such metamaterials may have wide range of nontrivial shapes of 
isofrequency contours including lemniscate, diamond and multiply connected curves with connectiv¬ 
ity number reaching five. The possibility of the dispersion engineering paves a way to more flexible 
manipulation of electromagnetic waves. Employing first-principle considerations we prove that uni¬ 
axially polarizable media should be described in terms of the nonlocal permittivity tensor which by 
no means can be reduced to local permittivity and permeability even in the long-wavelength limit. 
We introduce an alternative set of local material parameters including quadrupole susceptibility 
capable to capture all of the second-order spatial dispersion effects. 
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There exists a class of uniaxially polarizable media 
(UPM) that includes metamaterials based on coupled 
rod pairs, dual wires, etc. M exhibiting such inter¬ 
esting properties as artificial magnetism and backward 
waves |jl5[. The basic idea of coupled rod pair meta¬ 
material design is rather simple. The incident electro¬ 
magnetic wave excites currents in the coupled rods. The 
excitation of the symmetric mode when currents have 
the same magnitude and direction is usually attributed 
to the electric dipole resonance of the unit cell (Fig.[lji). 
The excitation of the antisymmetric mode when currents 
have the same magnitude but the opposite directions is 
attributed to the magnetic dipole resonance of the unit 
cell (Fig. [T|o) . Thus, one may conclude that the struc¬ 
ture possesses resonant artificial dielectric and magnetic 
properties due to the presence of the symmetric and anti¬ 
symmetric modes of the unit cell, respectively. Effective 
medium model (EMM) for coupled rod pair metamateri¬ 
als and similar structures in terms of local permittivity 
f loc and permeability p} oc tensors was introduced in the 
theoretical works Refs. verified by numerical sim¬ 

ulations [9] and by experiments [Tl fT()lIbl based on the 
Nicolson-Ross-Weir extraction procedure pd] Hi- 

In this Communication, we demonstrate that the exact 
description of the entire class of uniaxially polarizable 
media is incompatible with the local effective medium 
model, introduce an alternative set of local material pa¬ 
rameters and reveal the unusual dispersion properties of 
UPM. The isofrequency contours that are lemniscate- 
like, diamond-shaped and multiply connected (with con- 
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FIG. 1: (Color online) The sketch of (a) symmetric and (b) 
antisymmetric modes in coupled rod pair structure. Arrows 
show the direction of the currents. 


nectivity number equal to five) appear if d 2 e zz /dk 2 > 
0, £ zz (oj, 0) > 0 and 0 < d 2 e zz /dk 2 < 2 c 2 /uj 2 or 
d 2 e zz /dk 2 < 0, respectively. 

The nonlocal permittivity tensor fl6UT7l of UPM con¬ 
sisting of reciprocal elements with inversion symmetry 
can be reduced to22 


G 0 

f(w, k) = | 0 1 

t0 0 e zz 



(1) 


where the 2 axis is chosen in such way that P = P z z. 
Within the frame of EMM, the effective nonlocal permit¬ 
tivity tensor can be written in terms of e loc and ft} oc as 
follows EH: 


£eff (w, k) = e loc (w) + 2 k x (V oc (w)) 1 - / ) x k , (2) 


where q = uj/c, c is the speed of light, and the bian- 
isotropic terms are excluded since they vanish due to 
the structure inversion symmetry. Comparing Eq. m 
and Eq. |[2|) it is easy to show that uniaxially polariz¬ 
able medium is described by EMM in the only case when 

= loc = /,30 
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FIG. 2: (Color online) Illustration of the difference in mag¬ 
netic response of SRR and coupled rod pair. 
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In the other words, local permeability different from 
unity cannot be attributed to uniaxially polarizable me¬ 
dia. The qualitative explanation of the obtained result 
is as follows. We consider the excitation of the struc¬ 
ture (Fig. [2]) by the waves of two types that have the 
same direction of magnetic field but different directions 
of electric field and wave vector. In the case of split-ring 
resonator (SRR), the response to these two excitations 
is determined mainly by the direction of magnetic field 
and therefore is almost the same. On the other hand, 
coupled rod pair responds differently to these excitations 
as the alignment of the electric field vector with respect 
to the pair axis is crucial. This simple reasoning suggests 
that the magnetic response of coupled rod pair structure 
cannot be described properly by the local permeability 
tensor. 

Thus, it is important to introduce an alternative set 
of local parameters giving the UPM consistent descrip¬ 
tion. From now on we assume that for the UPM descrip¬ 
tion it is sufficient to consider only spatial dispersion of 
the second order and that higher order nonlocal effects 
are negligible. Under this approximation we may expand 
£zz{w,k) in Taylor series with respect to k. Due to the 
presumed inversion symmetry the first-order terms of the 
expansion vanish and it acquires the form: 

£zz(w, k) = £ loc (w) + Xij(u) hkj , (3) 

where e loc = £ 22 (w, 0). The quadrupole susceptibility 

I o2 

matrix Xij (<*>) = 5 d^dk depends only on the frequency 
and is symmetric with respect to the indices i and j. 
Thus, the quadrupole susceptibility \ yields an alterna¬ 
tive set of local material parameters, the number of inde¬ 
pendent components being equal to six due to the matrix 
X symmetry (the same as in the symmetric permeability 
tensor). Using Eq. ©, it is easy to obtain the relation for 
the macroscopic fields in space-time domain applicable to 
the general UPM: 

Dx - E X , Dy - Ey , 

D z = £ loc E Z ~Y2 Xij didj E z . (4) 

The symmetries of the structure reduce the number of 
the matrix x independent components fl9| . The nonzero 
components of quadrupole susceptibility for structures 
of different symmetries are summarized in the Table [T] 
For the sake of simplicity we consider the UPM with the 
symmetry C 2 , i.e. the structure invariant under rotation 
at the angle n around the anisotropy axis z. For example, 
this is the case for the coupled rod pair structure. The 
quadrupole susceptibility is given by the matrix 


X = 


Xxx 

Xxy 

V 0 


Xxy 0 \ 
Xyy 0 
0 Xzz) 


(5) 


Symmetry 

Nonzero components of x 

Independent 

components 

No 

Xxx, \yy 5 Xzz 5 Xxy = Xyx-, 
Xxz = Xzx 5 Xzy = Xyz 

6 

c 2 

Xxx1 Xyyi Xzz1 Xxy = Xyx 

4 

C 3 or C 4 

Xxx — Xyy> Xzz 

2 


TABLE I: The nonzero components of quadrupole suscepti¬ 
bility for the structures of different symmetries 


where the components depend on the geometry of the 
particular material and on the frequency. Now we con¬ 
sider the dispersion equation for TM waves in UPM 

k x + ky = ( q 2 — kl) £ zz (cc, k) (6) 

with £ 22 (w, k) given by Eqs. ©, ([5]). To illustrate phys¬ 
ical effects possible in UPM, we study the propagation 
in the plane Oyz when k x = 0, and in this case Eq. © 
reduces to 

k l [ X - (<? 2 - k l) Xyy] = (q 2 - h\) (£ loc + Xzzkl) . (7) 

The equation of the similar structure is valid for 
the plane of propagation k x = k y = kj_/y/2. In 
such case, k y should be replaced by k ±, and Xyy by 
1/2 ( Xxx + Xyy + 2 Xxy) in Eq. ©. 

In conventional local media with anisotropic uniaxial 
permittivity and permeability tensors, one can obtain 
elliptic and hyperbolic isofrequency contours, i.e. the 
second-order curves [2(J. Equation © yields that the 
isofrequency contours in UPM are the fourth order curves 
with the connectivity number more than two. This paves 
a way to more flexible electromagnetic waves manipula¬ 
tion using the nonlocal effects. Solving Eq. © we can 



FIG. 3: (Color online) The typical isofrequency contours for 
uniaxially polarizable medium with y 22 < 0. The plane k x = 
0 is examined. 
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FIG. 4: (Color online) The typical isofrequency contours for 
uniaxially polarizable medium with \ zz > 0. The plane k x ~ 
0 is examined. 



FIG. 5: (Color online) The geometry of the extraction exper¬ 
iment for uniaxially polarizable medium. TM polarization of 
the incident wave is explored. Coupled rods are parallel to 
the interface. 


visualise the dispersion properties of UPM. The system 
of isofrequency contours for various values of e , Xzz 
and Xyy is plotted in Figs. [3] and |U As one can see, 
UPM’s exhibit interesting and unusual dispersion prop¬ 
erties. Among the curves presented in Figs. 0 [4] there 
are both isofrequency contours that have already been 
observed for particular systems and the isofrequency con¬ 
tours that to our knowledge have not been observed so 
far for any particular structure. Namely, the isofrequency 
contour Fig. [3] for e loc > 0 and 0 < Xyy < 1 /q 2 is the 
same as one in the low-e mixed dispersion regime in the 
structure composed of the uniaxial electric scatterers [ 2 l| . 
Curves Fig.[3]for Xyy < 0 have a similarity with the isofre¬ 
quency contours obtained for the structure composed of 
dielectric rods J22|. A lemniscate-like isofrequency con¬ 


tour arises in the regime Xzz > 0 (Fig. [4]), e loc ss 0 and 
0 < Xyy < l/<? 2 - The structures with such dispersion 
properties can be useful for obtaining superradiance or 
for phase matching (if nonlinear elements are present). A 
diamond-shaped isofrequency contour appears if Xzz > 0 , 
e(tu, 0) > 0 and 0 < Xyy < 1 /g 2 ; such dispersion proper¬ 
ties are useful in Purcell factor engineering [23]. 

It is also important that the dispersion equation Eq. 0 
is linear with respect to k 2 and biquadratic with respect 
to k z . Therefore, if a surface of the structure is normal 
to y axis (Fig [5]), a single transmitted beam will be ob¬ 
served and the spatial dispersion effects are hidden to 
some extent. However, if the surface of the structure 
is perpendicular to the anisotropy axis z, two transmit¬ 
ted TM modes will be excited which serves as a direct 
manifestation of the nonlocal effects. 

To determine the type of the dispersion regime for the 
particular UPM one has to infer the quadrupole suscep¬ 
tibility matrix x of the structure. This can be fulfilled 
employing the technique similar to Nicolson-Ross-Weir 
method |l3fl5| . In the case k z = 0 (Fig. [5]) Eq. J7J) re¬ 
duces to the following one: 



£ loc 

1 - 9 2 Xyy ' 


( 8 ) 


Comparing this result to the standard formula ky = 
q 2 £ loc y} oc we conclude that it is possible to introduce 
local “permeability” for this particular direction of prop¬ 
agation: 


A* 


loc 


1 

1 - 9 2 Xyy ' 


(9) 


If the function on the right-hand side of Eq. (0 is decreas¬ 
ing with the frequency, the wave described by Eq. 0 is 
backward. It is this backward wave that was observed in 
coupled rod pair metamaterials [Ej]. Moreover, assuming 
the sectional continuity of the fields at the UPM sur¬ 
face and integrating Maxwell equations over a pillbox, it 
can be easily demonstrated that the boundary conditions 
for the TM wave incident at the UPM surface from vac¬ 
uum are equivalent to those for the boundary between 
magnetic 0 and vacuum in the special case of normal 
incidence ( k x = k z = 0 in the geometry Fig. 0: 


E\ z — E 2z , 

B lx = (/r loc ) _1 B 2x . 

that corresponds to the boundary conditions for the 
conventional magnetic with the local parameters e loc 
and /i loc , and the reflection coefficients determined in 
the case of normal incidence will coincide with those 
for the magnetic Eq. 0. Thus, e loc and y} oc can 
be extracted by the standard Nicolson-Ross-Weir tech¬ 
nique HEHm, and Xyy can be calculated via y} oc as 
Xyy = (l — (/i loc ) -1 ) /q 2 . The analogous approach is 
valid for Xxx extraction; for this purpose the plane of the 
pairs shold be parallel to the boundary of the structure. 
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However, in the case of oblique incidence with k z ^ 0 
EMM is no longer valid. From the measured reflection 
at oblique incidence it is possible to infer \zz- Investi¬ 
gation of the reflection of the TM wave incident in the 
plane different from Oyz and Oxy (namely, the plane 
k x = k y ) will allow one to extract \xy component of the 
structure quadrupole susceptibility. The detailed discus¬ 
sion of UPM quadrupole susceptibility extraction will be 
published elsewhere. 

The reasoning similar to the performed above reveals 
that for the structures that have the effective nonlocal 
permittivity tensor of the form 


e(cj, k) = 



( 11 ) 


the only possible nontrivial component of the magnetic 
permeability tensor is y zz component, i.e. p = xx + 
yy + y zz zz. It is interesting to notice that the fishnet 
structure Q made from reciprocal material comprises the 
class of materials described by Eq. (EH). Indeed, e zx = 
E zy = 0 and s zz = 1 for the fishnet structure that can be 
polarized in the plane Oxy. Making use of the symmetry 
of the effective permittivity tensor, we obtain also that 
e xz = £yz = 0. Therefore, quadrupole susceptibilities 
should be also used for the description of the structures 
Eq. En¬ 
in this work, we reveal the intriguing dispersion prop¬ 
erties of uniaxially polarizable media, namely, nontriv¬ 
ial shapes of isofrequency contours including lemniscate, 


diamond and multiply connected curves with connectiv¬ 
ity number reaching five. All these unusual properties 
are beyond effective medium model. Problems arising 
in the effective medium approach were illustrated in a 
number of works f2dl [2?| . In the present Communication, 
we demonstrate incompleteness of local effective medium 
model for the entire class of uniaxially polarizable media 
proving from the first principles that these structures do 
not possess magnetic properties (/t loc = J) though they 
do exhibit second-order spatial dispersion effects. We 
also propose an alternative set of local material param¬ 
eters that are capable to capture all of the second-order 
nonlocal effects. This set of parameters allows one to 
demonstrate how the backward waves arise in uniaxially 
polarizable media due to the spatial dispersion effects 
and provide a correct interpretation to the experiments 
extracting material parameters of uniaxially polarizable 
media. 
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We assume here that e(oj,—k) = e(oj,k), i.e. structure un¬ 
der the study possesses inversion symmetry. Under this 
assumption and employing the symmetry of the kinetic co¬ 
efficients 0111 , it can be shown that £ik{ca, k) = £ki( w, k), 
i.e. the effective nonlocal permittivity tensor is symmetric. 
Therefore, £ zx = £zy = 0. As a result, the only nontrivial 
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30 Indeed, comparing Eqs. ©> © we obtain k x A x k = 
A (fc) zz where z denotes a unit vector along the z axis, A = 
^(/l loc (a;)) _1 — I^j /q 2 , and the term A (k) = £ zz (u,k) — 
£ zz (oj, 0) is the spatial dispersion corrections to £ zz . This 
matrix equation yields eight scalar equations with the zero 
right hand side. Taking into account that the wave vector 
components k x , k y and k z are independent, one obtains 
a homogeneous system of nine linear equations with the 
nonzero determinant. Thus, we infer that the only possi¬ 
ble solution for the matrix I is a zero matrix and, as a 
consequence, p}° c = I. 




